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Abstract:  Mathematical  expressions  have  been  estab¬ 
lished  for  estimating  the  bulk  salinity  of  Arctic  and 
Antarctic  sea  ice  vs.  ice  floe  thickness.  The  ice  salinity 
vs.  thickness  relationships  are  based  on  data  for  over 
400  sea  ice  cores  compiled  from  numerous  sources. 
The  results  show  that  the  bulk  salinity  of  first-year  sea 


ice  decreases  in  an  exponential  trend  with  ice  sheet 
thickness.  A  similar  trend  reoccurs  as  the  winter  ice 
passes  through  the  melt  season.  The  expression  for  the 
bulk  salinity  Sg  in  %o  for  first-year  sea  ice  from  1 0  to 
200  cm  thick  is  =  4.606  -i-  91 .603/rp,  where  is 
the  ice  floe  thickness  in  centimeters. 


Cover:  Logging  of  a  sea  ice  core  in  the  field.  After  extraction  from  the  ice  sheet,  the  ice  core  is  quickly  removed 
from  the  ice  core  barrel  and  placed  on  a  cutting  board,  not  on  the  snow,  which  can  draw  brine  from  the  ice.  The 
core  length  is  measured.  Small  holes  are  then  drilled  about  2-3  cm  deep  into  the  core  side  at  20-cm  increments 
and  a  thermistor  is  inserted  in  each  hole  to  determine  the  in-situ  temperature  of  the  ice  siieet.  Closer  spaced 
measurements  may  be  required  near  the  ice  sheet  surface,  top  of  the  ice  core,  to  better  define  the  temperature 
gradient  in  the  region.  The  core  is  next  cut  into  lengths  using  the  cut-off  saw,  which  is  mounted  at  the  end  of  the 
cutting  board,  and  the  axial  length  of  the  ice  cylinder  then  is  determined  (in  this  case,  using  a  dial  gauge 
measurement  jig  located  beside  the  notebook).  The  ice  sample  is  next  weighed  on  an  electronic  balance  and 
then  placed  in  a  sealed  plastic  bag  or  container  for  later  melting  and  salinity  measurements.  Air  and  ice  sheet 
temperature,  ice  porosity  and  processing  time  are  factors  that  effect  the  loss  of  brine  from  the  ice  core  during 
logging.  Upward  flow  of  seawater  through  the  permeable  ice  remaihing  at  the  bottom  of  a  borehole  will  flush 
brine  from  the  ice.  When  this  ice  is  finally  removed  by  coring,  its  brine  concentration  may  have  chohged 
appreciably  and  the  subsequent  salinity  measurement  would  be  in  error. 
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INTRODUCTION 

When  seawater  freezes,  small,  flat  platelets  of 
fresh  ice  are  first  to  form.  As  this  ice  mass  grows, 
the  platelets  commingle  and  bond  one  to  another 
to  form  a  compliant,  highly  saline  slush  called 
grease  ice.  With  sustained  freezing  temperatures, 
the  platelet  mass  thickens  and  grows  together, 
forming  an  ever-stiffening  ice  cover.  During  this 
ice  growth  process  most  of  the  seawater  is  dis¬ 
placed,  but  between  the  platelets  of  pure  ice  a 
seawater  concentrate  called  brine  becomes  en¬ 
trapped.  This  entrapment  causes  the  ice  to  be 
salty:  the  salinity  of  the  surface  layer  of  the  new 
ice  can  be  greater  than  40  %o  (Martin  1979,  Kasai 
and  Ono  1984,  Ono  and  Kasai  1985). 

Many  factors  affect  the  thickness  of  the 
platelet  or  frazil  ice  layer,  including  growth 
rate,  sea  state,  ice  pack  motion,  and  under-ice 
turbulence.  However,  under  relatively  quies¬ 
cent  conditions  and  after  the  frazil  ice  layer 
has  "stabilized,"  the  stage  is  set  for  the  growth 
of  congelation  ice.  At  first,  the  growing  ice 
platelets,  <0.5  mm  thick,  are  randomly  ori¬ 
ented,  but  in  time  they  become  more  and  more 
ordered,  grow  wider,  and  stack  one  against 
another.  Tliey  freeze  together  in  layered  groups 
to  form  individual  ice  crystals.  Due  to  thermo¬ 
dynamic  factors,  preferred  platelet  growth  in 
seawater  is  vertical  (Weeks  and  Ackley  1989). 
Therefore,  over  a  short  growth  span  the  ice 
crystals  become  "vertically"  aligned,  with  the 
parallel  platelet  array  having  a  predominantly 
horizontal  c-axes  alignment.  The  ensuing 
growth  structure  is  called  congelation  ice.  The 
number  and  size  of  the  vertical  ice  platelets  in 
any  crystal  vary  with  growth  rate,  as  does  the 
volume  of  brine  trapped  between  the  fresh-ice 
platelets  (Cox  and  Weeks  1988,  Weeks  and 
Ackley  1989). 


Brine  entrapment  originates  in  the  root  area 
(Fig.  1)  of  the  1-  to  5-cm-thick  mushy  region  called 
the  SK  or  skeleton  layer  (Harrison  and  Tiller  1963) 
or  crystallization  front  (Cherepanov  1975).  In  this 
growth  zone  the  individual  ice  platelets  extend 
into  the  seawater,  like  fingers  on  a  hand,  and  are 
typically  less  than  0.25  mm  thick  and  more  than  1 
cm  wide.  As  each  dendritic  platelet  grows  thicker 
at  its  root,  it  freezes  to  its  neighbor.  A  high  growth 
rate  at  this  stage  (>2  cm/ day)  results  in  platelets 
less  than  0.5  mm  thick,  but  at  low  growth  rates 
(<0.5  cm/day)  the  platelets  can  grow  to  be  1.6 
mm  thick,  with  a  related  increase  in  crystal  size 
(Weeks  and  Hamilton  1962a,b;  Page  1966;  Nakawo 
and  Sinha  1984).  During  this  crystal  formation 


a  <  b  <  c 

a  ~  0.1  to  0.3  mm;  b  ~  1  to  5  x  a;  c  >  5  x  a 
d  ~  0.25  to  1 .25  mm  {avg  0.7) 

Frozen  Interlace 
I’-'  I  Seawater  Interface 

Figure  1.  Sea-ice  crystal  structure. 


Figure  2.  Sea-ice  core  sections  (from  Morey  et  al  1984). 


process  the  seawater  at  the  ice  growth  interface  is 
replaced  by  supersaturated  brine  that  slows  the 
rate  of  solidification.  As  with  the  frazil/ granular 
surface  ice  layer,  some  brine  and  gas  become 
trapped  between  the  platelets,  as  depicted  in  Fig¬ 
ure  1.  As  the  ice  sheet  continues  to  grow  and 
becomes  colder,  the  film  of  brine  between  the  ice 
platelets  reduces  in  volume  due  to  expulsion  and 
drainage  processes  as  well  as  the  freezing  out  of 
fresh  water  from  the  brine  to  the  surrounding  ice 
surface.  The  now-isolated  brine  inclusions  are 
called  brine  pockets.  They  are  frequently  described 
as  ellipsoids,  but  in  fact  can  have  rather  complex 
geometries,  depending  on  ice  growth  rate  and 
temperature. 

Brine  expulsion  and  drainage  occur  along  a 
flow  route  that  resembles  the  trunk  and  branches 
of  a  tree,  called  a  drainage  tube  and  channels, 
respectively.  A  vertical  drainage  tube  and  chan¬ 
nel  system  in  Beaufort  Sea  congelation  ice  is 
shown  in  Figure  2.  From  1.5-m-square  ice  blocks 
removed  from  1-  and  1.8-m-thick  ice  sheets  in  the 
Beaufort  Sea,  I  found  the  drainage  tubes  at  the 
bottom  of  the  ice  to  be  spaced  5  to  17  cm  apart, 
but  they  average  10  cm  apart.  These  spacings  are 
in  agreement  with  the  findings  of  Martin  (1979), 
Wakatsuchi  (1983),  and  other  investigators. 
Wakatsuchi  also  reported  that  the  spatial  density 
of  brine  tubes  depends  on  the  growth  rate.  The 
faster  the  growth  rate,  the  more  drainage  tubes 
per  unit  area.  The  brine  channel  inclination  varies 
from  40°  to  54°,  but  averages  45°  (Lake  and  Lewis 
1970,  Criminale  and  LeLong  1984), 

Brine  drainage  is  high  during  initial  ice  for¬ 
mation,  but  as  the  ice  sheet  thickens  and  becomes 


colder,  brine  flow  decreases  and  much  of  the  drain¬ 
age  system  is  transformed  into  a  series  of  brine 
and  gas  cells  separated  by  necks  of  ice  (Doronin 
and  Kheisin  1975).  The  salinity  at  any  given  depth 
now  decreases  slowly  during  the  remainder  of 
the  growth  season.  Drainage  systems  form  again 
during  the  melt  season  (Nakawo  and  Sinha  1981). 
Near-freezing  temperatures  and  solar  radiation 
heating  cause  melting,  "corrosion  pockets"  form 
in  the  ice  surface,  and  internal  melting  occurs  at 
crystal  boundaries,  old  brine  drainage  routes, 
flaws,  etc.  (Bermington  1967,  Cox  and  Schultz 
1980).  The  corrosion  pockets  contain  liquid  with  a 
high  salt  concentration  originating,  for  example, 
from  the  salt  residue  of  the  >50  %o  salinity  brine 
that  was  ejected  to  the  surface  during  initial  ice 
sheet  formation  (Burke  1940,  Zubov  1945,  Martin 
1979,  Drinkwater  and  Crocker  1988,  Perovich  and 
Richter-Menge  1994),  from  the  typically  high-sa¬ 
linity  frazil  or  granular  ice  surface  layer  (Martin 
1979)  or  from  brine  concentrated  in  an  ice  surface 
layer  formed  as  a  result  of  ice  sheet  flooding.  The 
enriched  brine  then  melts  into  the  ice  sheet,  re¬ 
opening  or  enlarging  the  old  drainage  systems  as 
the  brine  wave  descends  to  the  sea.  The  cold, 
high-salinity  (60  to  70  %o)  drainage  (Lewis  and 
Milne  1977)  is  followed  by  increasingly  fresher 
flow  resulting  from  the  pure-ice  melt  or  is  re¬ 
placed  by  an  upward  flow  of  less  dense  seawater 
(Edie  and  Martin  1975).  Depending  on  the  brine 
temperature  and  volume  flux  and  seawater  cur¬ 
rent  velocity,  on  exiting  the  ice  the  brine  drainage 
can  freeze  the  seawater  around  the  periphery  of 
the  descending  brine  plume,  forming  long,  deli¬ 
cate,  thin-walled  hollow  ice  stalactites  (Fig.  3)  that 
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Figure  3.  Over  1-m-long  stalactite 
that  formed  as  a  result  of  cold  brine 
drainage  into  the  sea  at  the  bottom 
of  an  icefloe. 


on  occasion  can  extend  up  to  6  m  below  the  bot¬ 
tom  of  the  sea  ice  (Page  1970,  Dayton  and  Martin 
1971,  Martin  1974,  Lewis  and  Milne  1977, 
Grishchenko  1988). 

As  reported  in  Kovacs  and  Morey  (1980),  the 
growth  of  congelation  sea  ice  in  cahn  water  re¬ 
sults  in  a  structure  in  which  the  c-axes  of  the 
crystals  are  randomly  oriented  in  the  horizontal 
plane.  This  is  depicted  in  Figure  4.  However,  if  a 
current  exists  under  the  ice  sheet,  selective  ice 


C-axis 

Orientation 


Parallel  Fresh  Ice 
Platelets 


Brine  and  Gas 
Pockets 


Figure  4.  Three  columnar  sea- 
ice  crystals  with  unaligned  c- 
axes.  The  crystal  boundaries 
are  highly  irregular. 


platelet  growth  and  therefore  crystal  fabric  orien¬ 
tation  occurs,  in  which  the  horizontal  c-axes  of 
the  ice  crystals  become  aligned  with  the  current 
(Kovacs  and  Morey  1978,  Weeks  and  Gow  1978, 
1979;  Langhome  and  Robinson  1986).  The  signifi¬ 
cance  of  this  alignment  is  that  it  renders  the  ice 
sheet  horizontally  anisotropic,  and  this  in  turn 
affects  the  electromagnetic  (Campbell  and  Orange 
1974;  Kovacs  and  Morey  1978, 1979;  Morey  et  al. 
1984)  and  mechanical  (Anderson  1963,  Payton 
1966,  Wang  1979,  Timco  et  al.  1991,  Kovacs  1993) 
properties  of  the  ice.  Many  other  investigators 
have  observed  preferred  c-axis  alignment  in  sea 
ice.  Perhaps  most  notable  is  Cherepanov  (1975), 
who  made  wide-area  sea-ice  observations  in  the 
seas  north  of  Russia,  and,  more  recently,  Gow  et 
al.  (1982),  who  verified  that  similar  sea-ice  c-axis 
alignment  occurs  in  Antarctic  congelation  sea  ice. 

While  the  above  description  of  first-year  sea- 
ice  growth  does  not  encompass  the  many  struc¬ 
ture  variations  that  can  occur  due  to  such  factors 
as  snow  loading  and  flooding,  thick  frazil-ice 
components  (e.g.,  Cherepanov  and  Kozlovskiy 
1973a,b;  Martin  1981;  Gow  et  al.  1987),  and  large 
interlaced  ice  platelet  layers,  as  found,  for  ex¬ 
ample,  in  Antarctic  sea  ice  (e.g.,  Cherepanov  and 
Kozlovskiy  1973a,  Dieckmann  et  al.  1986,  Eicken 
and  Lange  1989,  Jeffries  and  Weeks  1993),  it  does 
reveal  that  sea  ice  is  a  multicomponent  material. 
It  is  composed  of  fresh  ice,  brine,  and  gas  inclu¬ 
sions  and  solid  salt  crystals  (Sinha  1977,  Marion 
and  Grant  1994).  The  volume  of  fresh  ice  is  by  far 
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the  largest  fraction,  typically  in  excess  of  95%.  Sea 
ice  is  classified  by  age  (first-year,  second-year, 
and  multiyear)  and  by  morphology  (flat,  ridge, 
rafted,  etc.).  Variations  in  growth,  melt,  and  de¬ 
formation  processes  result  in  ice  formations  of 
complex  shape  and  structure  and  variable  brine 
and  gas  content.  The  sea  ice  may  also  contain 
foreign  matter  such  as  sediment  (e.g.,  Barnes  et  al. 
1984,  Kempema  et  al.  1989)  and  microscopic  algae 
and  diatoms  (e.g.,  Sutherland  1852;  Wright  and 
Priestley  1922;  Usachev  1949;  Buinitskii  1965, 1974; 
and  many  others  more  recently). 

ICE  FLOE  BRINE  VOLUME 
AND  SALINITY 

A  number  of  factors  affect  the  volume  of  brine 
that  is  initially  incorporated  into  the  sea  ice.  Fore¬ 
most  is  the  salinity  of  the  seawater,  the  ice  growth 
rate  and  structure.  Under  extreme  cold  condi¬ 
tions,  accelerated  growth  occurs,  and  more  brine 
tends  to  be  entrapped  in  the  ice  at  the  bottom  of 
the  ice  sheet  (e.g.,  Malmgren  and  Sverdrup  1927, 
Burke  1940,  Zubov  1945,  Weeks  and  Lee  1958, 
Nakawo  and  Sinha  1981)  and  for  congelation  ice 
the  ice  platelets  that  form  each  ice  crystal  are 
narrower  (e.g..  Weeks  and  Hamilton  1962a,b;  Page 
1966;  Nakawo  and  Sinha  1984).  As  the  ice  sheet 
thickens,  the  rate  of  growth  slows  due  to  the  de¬ 
creasing  rate  at  which  heat  is  transferred  from  the 
bottom  of  the  ice  to  the  atmosphere  (Makshtas 
1984,  Maykut  1986).  This  heat  exchange  is  further 
reduced  with  increasing  snow  cover.  The  net  ef¬ 
fect  is  that  less  brine  entrapment  occurs  as  the  ice 
increases  in  thickness.  Over  time,  very  gradual 
brine  migration  occurs  toward  the  warmer  sea¬ 
water  (Whitman  1926,  Zubov  1945,  Kingery  and 
Goodnow  1963,  Bennington  1967,  Untersteiner 
1968,  Nakawo  and  Sinha  1981).  An  explanation 
for  this  continued  brine  migration  is  the  dissolu¬ 
tion  of  ice  by  the  brine  at  the  warmer  bottom  end 
of  a  brine  pocket  and  a  refreezing  of  water  at  the 
colder  upper  end  of  the  pocket.  This  migration 
and  associated  brine  drainage  gradually  reduce 
the  bulk  salinity  of  the  ice  from  an  initial  value  of 
as  much  as  25  %o  for  ice  5  cm  thick  to  about  5%o 
for  first-year  ice  about  2  m  thick. 

Accelerated  brine  drainage  can  occur  under 
warming  conditions.  The  enriched  brine  from  the 
top  of  the  ice  sheet  can  move  downward  like  a 
wave  slowly  descending  through  the  ice  sheet. 
During  the  descent  the  general  form  of  the  wave¬ 
let  is  preserved,  but  the  amplitude  (that  is,  the 


Figure  5.  Major  17 -cm  brine  wave  mi¬ 
gration  event  observed  in  Antarctic  fast 
sea  ice  between  14  May  and  5  June 
(after  Doronin  and  Kheisin  1975). 

brine  volume)  gradually  decreases.  This  is  due  to 
brine  pockets  intercepting  the  drainage  network 
and  draining  to  the  sea  in  advance  of  the  follow¬ 
ing  brine  wave.  This  drainage  phenomenon,  il¬ 
lustrated  in  terms  of  a  very  smooth  ice  sheet  sa¬ 
linity  distribution  profile,  is  shown  in  Figure  5. 
The  salinity  (brine)  migration  wave  shown  is  very 
fast,  about  0.8  cm/ day.  Migration  rates  one  to  two 
orders  of  magnitude  lower  may  be  more  typical 
for  winter  sea  ice,  based  on  the  velocity  assess¬ 
ments  presented  by  Weeks  and  Ackley  (1989),  who 
give  an  in-depth  review  of  sea-ice  salinity  and 
brine  entrapment  and  migration  processes.  In  ad¬ 
dition,  it  is  difficult,  at  best,  to  see  this  movement 
in  the  similar  data  of  Nakawo  and  Sinha  (1981), 
and  a  rate  of  1  mm/day  may  be  estimated  from 
Shapiro  and  Weeks  (1993,  Figure  9). 

A  widely  cited  study  on  the  bulk  salinity  varia¬ 
tion  in  sea  ice  is  that  of  Cox  and  Weeks  (1974).  In 
this  study  average  or  bulk  salinity  data  for  winter 
Beaufort  Sea  ice  was  collected  from  a  number  of 
sources.  Their  plot  of  the  salinity  vs.  ice  thickness 
data  clearly  shows  a  decrease  in  bulk  salinity  with 
increasing  ice  floe  thickness  (Fig.  6).  The  available 
data  seemed  to  justify  representing  the  bulk  sa¬ 
linity  vs.  ice  floe  thickness  trend  by  the  two  linear 
curves  shown  in  Figure  6.  However,  many  natu¬ 
ral  phenomena  tend  to  follow  exponential-like 
distributions.  Sea-ice  bulk  salinity  vs.  floe  thick¬ 
ness  appears  to  be  one  of  them. 
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Figure  6.  Sea-ice  bulk  salinity  vs.  floe  thickness  as  reported  by  Cox  and  Weeks  (1974). 


Figure  7.  Sea-ice  bulk  salinity  vs.  ice  sheet  thickness  and  growth 
time  (after  Weeks  and  Lee  1958). 


In  the  mid-1950s.  Weeks  and  Lee  (1958)  pub¬ 
lished  their  observations  on  the  bulk  salinity  of 
an  ice  sheet  vs.  ice  thickness  and  growth  time 
(Fig.  7).  The  abrupt  step  in  the  curve  through  the 
data  occurred  as  a  result  of  a  flooding  event  dur¬ 
ing  which  seawater  infiltrated  the  snow  cover. 
Had  this  not  occurred,  one  could  visualize  that 
the  trend  of  the  curve  before  the  step  would  have 
continued  to  decrease  to  a  salinity  of  about  5%o 
by  May.  In  any  event,  the  data  clearly  reveal  a 


rapid  decrease  in  bulk  salinity  as  the  ice 
sheet  grew  to  about  0.5  m  thick  and  a 
gradual  decrease  in  bulk  salinity  there¬ 
after. 

Other  early  reports  also  confirm  that 
this  bulk  salinity  vs.  floe  thickness  trend 
applies  to  Antarctic  sea  ice.  Fedotov 
(1973)  measured  the  bulk  salinity  of  the 
fast  ice  in  Alasheev  Bight  during  the 
growth  period.  His  results  are  presented 
in  Figure  8.  The  bulk  salinity  was  found 
to  be  about  25  %o  in  the  5-cm-thick  ice 
but  dropped  rapidly  to  about  7.5  %o  when 
the  ice  was  0.5  m  thick.  At  the  end  of  the 
measurement  period  (about  170  days), 
the  ice  was  150  cm  thick  and  the  bulk 
salinity  had  decreased  to  about  5%o. 
Doronin  and  Kheisin  (1975)  presented 
results  for  the  bulk  salinity  of  Antarctic 
fast  sea  ice  during  the  growth  and  melt  periods 
(Fig.  9).  This  figure  again  shows  a  rapid  decrease 
in  the  bulk  salinity  during  the  first  growth  month, 
followed  by  a  gradual  decrease  until  the  onset  of 
the  melt  season  in  October.  At  this  time  another 
rapid  decrease  in  the  bulk  salinity  began,  which 
ended  in  December.  The  ice  sheet's  bulk  salinity 
had  by  then  decreased  to  about  2%o,  a  value  rep¬ 
resentative  of  Arctic  multiyear  sea  ice  of  similar 
thickness. 
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Figure  8.  Sea-ice  hulk  salinity  and  growth  rate  vs.  floe  thickness 
(after  Fedotov  1973). 


Growth  Time 


Figure  9.  Sea-ice  bulk  salinity  vs.  growth  time  (after  Doronin 
and  Kheisin  1975). 


ESTIMATING 
ICE  FLOE  SALINITY 

Ryvlin  (1974)  was  one  of  the  first  to 
express  the  bulk  salinity  Sg  vs.  ice  floe 
thickness  Tp  trend  empirically.  His  for¬ 
mula  takes  into  consideration  the  salinity 
of  the  seawater  and  the  ice  sheet 
growth  rate  as  follows: 

SB=S^(l-SR)e-5(TFr +Sr  (1) 

where  Sr  =  Sgp/S^  =  salinity  ratio 

Spp  =  final  bulk  salinity  at  end  of 
growth  season 

a  =  growth  rate  coefficient,  re¬ 
ported  to  vary  from  0.35  to 
0.5,  e.g.,  at  Gr  >  4  cm/ day, 
a  =  0.35,  and  at  Gr  <  0.5 
cm/day,  a  =  0.60. 

Ryvlin  suggests  that  where  Gr  is  not 
known,  one  may  assume  the  value  of  0.5 
for  a.  He  also  found  empirically  that  Sr  is  ~ 
0.13.  Using  a  value  of  0.5  and  34  %o  for  a 
and  S^  respectively,  and  an  unknown  value 
for  Srp,  Ryvlin  compared  eq  1  with  his  and 
Weeks  and  Anderson's  (1958)  first-year  sea 
ice  Sr  vs.  Tp  field  data,  as  shown  in  Figure 
10.  In  this  comparison,  the  agreement  be¬ 
tween  the  empirical  estimate  and  field  mea¬ 
surements  is  extremely  good  throughout 
the  entire  ice  floe  thickness  range.  The  ques¬ 
tion  is,  does  this  agreement  hold  for  other 
data  sets? 

To  evaluate  Ryvlin's  equation,  we  first 
compiled  unpublished  first-year  Beaufort 
Sea  ice  Sr  vs.  Tp  data  that  was  provided 
by  Perovich  (five  cores  for  Tp  <25  cm) 


The  decrease  in  bulk  salinity  vs.  time  shown  in  and  Kovacs  (30  cores),  along  with  published  data 


Figures  7  and  9  was  also  observed  during  lab-  from  Gow  and  Weeks  (two  cores,  1977),  Tucker  et 

oratory  studies  of  brine  entrapment  processes  al.  (five  cores,  1984),  and  Meese  (nine  cores,  1989) 

in  NaCl  ice  (Cox  and  Weeks  1975)  and  in  natu-  (Fig.  11).  The  salinity  of  the  Beaufort  Sea  near  the 

ral  sea  ice  by  Fukutomi  et  al.  (1951)  and  Blinov  Alaska  coast  is  31.5%o.  The  purpose  of  using  data 

(1965).  from  just  one  location  is  that  the  salinity  of  the 

These  reports  clearly  indicate  that  the  bulk  sa-  water  is  a  parameter  that  affects  ice  salinity  (Weeks 

linity  of  natural  and  laboratory-grown  saline  ice  and  Lofgren  1967,  Ryvlin  1974).  This  can  be  seen 

decreases  rapidly  with  increasing  thickness  or  time  in  Figure  12,  where  the  only  parameter  varied  in 

during  the  early  growth  phase,  followed  by  a  eq  1  was  the  seawater  salinity.  The  fbced  param- 

second  rapid  decrease  in  the  bulk  salinity  during  eters  used  in  eq  1  to  create  the  curves  in  Figure  12 

the  early  part  of  the  melt  season.  Because  of  the  were  0.5  for  a  and  0.13  for  Sr.  For  an  end-of- 

latter  decrease,  the  bulk  salinity  data  for  "warm"  growth-season  ice  thickness  of  200  cm,  Sr  is 

or  "old"  sea  ice  should  not  be  combined  and  ana-  4.44 %o  and  4.12%o  for  an  S^  of  34 %o  and  31.5 %o, 

lyzed  with  "cold"  first-year  sea  ice  data.  respectively.  For  this  example,  the  seawater  sa- 
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Figure  10.  Comparison  of  Ryvlin's  empirical  equation  for  esti¬ 
mating  sea-ice  bulk  salinity  vs.  floe  thickness,  with  field  data 
(after  Ryvlin  1974). 


Figure  11.  Beaufort  Sea  first-year  ice  bulk  salinity  vs.  floe  thick¬ 
ness.  The  lower  curve  represents  the  results  of  Ryvlin's  equa¬ 
tion  for  31.5%o  seawater. 


Unity  difference  caused  an  8%  difference  in  the 
bulk  sea-ice  salinity. 

The  regression  curve  through  the  Beaufort  Sea 
Sg  vs.  Tp  data  (Fig.  11)  is,  for  all  intents  and  pur¬ 
poses,  horizontal  after  an  ice  thickness  of  about 
150  cm.  Therefore,  Sgp  may  be  determined  from 
the  regression  curve  at  any  Tp  after  150  cm  with 
minor  effect  on  the  results  of  eq  1.  For  consis¬ 
tency,  SpF  at  Tp  =  200  cm  is  used  in  this  report. 

Running  below  the  regression  line,  and  most 
of  the  data,  in  Figure  11  is  the  curve  derived  using 
eq  1  and  an  of  31.5%o  (Fig.  12).  Even  though 
the  number  of  data  points  is  limited  (and  non¬ 
existent  for  a  Tp  of  25  to  50  cm)  it  is  apparent  that 


eq  1  does  not  fit  the  data.  This  is  surpris¬ 
ing,  given  the  good  agreement  shown  in 
Figure  10. 

At  200  cm  the  regression  curve  in  Fig¬ 
ure  11  gives  a  value  of  5.50%o  for  Sg.  Divid¬ 
ing  this  salinity  by  =  31.5  %o  gives  a 
value  for  Sr  of  0.175,  not  0.13  as  Ryvlin 
determined  from  his  data  analysis.  Using 
this  new  Sr  value  along  with  his  suggested 
growth  rate  coefficient  of  0.5  in  eq  1  brings 
the  Ryvlin  curve  into  good  agreement  with 
the  regression  curve  between  Tp  =  125  and 
200  cm.  However,  there  is  still  a  significant 
difference  between  eq  1  and  the  regression 
curve  where  Tp  is  less  than  about  50  cm.  To 
bring  the  two  curves  into  agreement,  coef¬ 
ficient  a  has  to  vary  with  Tp.  This  is  shown 
in  Figure  13,  where  the  lower  curve  repre¬ 
sents  eq  1  with  only  Sr  changed  from  0.13 
to  0.175  and  the  upper  curve  represents 
both  the  regression  line  through  the  data 
in  Figure  11  and  eq  1  with  Sr  =  0.175  and 
coefficient  a  varied  as  shown. 

To  verify  the  need  to  vary  Sr  and  in  eq 
1, 1  analyzed  the  Sg  vs.  Tp  data  collected  by 
Jeffries  (1994)  in  the  Amundsen  and 
Bellingshausen  Seas,  Antarctica,  in  late 
winter.  Unlike  the  predominantly  congela¬ 
tion  sea-ice  data  from  the  Arctic  (Beaufort 
Sea)  the  Antarctic  sea-ice  data  are  from  ice 
sheets  having  a  predominantly  frazil,  large 
interlaced  platelet  and  infiltrated  snow  ice 
structure.  Core  data  on  ice  up  to  200  cm 
thick  with  a  mean  temperature  below  -3°C 
were  selected.  The  latter  ensured  that  the 
ice  was  "cold"  and  not  melting.  Their  data 
from  74  ice  cores  are  shown  in  Figure  14. 
The  representative  regression  curve 
through  the  data  is  of  the  same  form  as  the 
one  through  the  Beaufort  Sea  ice  data  in 
Figure  11.  It  is  important  to  note  that  the  re¬ 
gression  curve  at  Tp  =  200  cm  gives  Sg  =  5.41 
vs.  5.50 %o  for  the  Beaufort  Sea  data.  Given  the 
fact  that  the  Antarctic  seawater  salinity  is  34  %o* 
vs.  31.5%o  for  the  southern  Beaufort  Sea,  one 
would  expect  these  Sg  values  to  be  reversed  for 
reasons  previously  discussed.  This  unexpected 
result  suggests  that  the  sea-ice  data  for  the 


M.  Jeffries  (1995),  Geophysical  Institute,  Univer¬ 
sity  of  Alaska,  Fairbanks,  and  S.  Jacobs  (1995), 
Lamont-Doherty  Geological  Observatory  of  Colum¬ 
bia  University,  Palisades,  New  York,  personal  com¬ 
munications. 
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Figure  12.  Sea-ice  bulk  salinity  vs.  floe  thickness 
as  calculated  with  Ryvlin's  equation  for  31.5  and 
34%o  seawater. 


Figure  13.  Ryvlin's  equation  for  calculating  sea- 
ice  bulk  salinity  vs.  floe  thickness  where  param¬ 
eter  a  is  constant  (lower  curve)  and  varied  (up¬ 
per  curve)  to  match  the  regression  curve  through 
the  Beaufort  Sea  ice  data  in  Figure  11. 


Arctic  (Fig.  11)  and  Antarctic  (Fig,  14)  are  not 
sufficient  to  fully  define  the  S^  vs.  Tp  trends  or 
that  the  effect  of  the  seawater  salinity  differ- 
ence  on  S5  is  masked  by  the  scatter  in  the  data 
and/or  ice  structure  differences.  Perhaps  as  im¬ 
portant  is  the  regression  curve  selected  to  repre¬ 
sent  the  data.  Statistically,  two  regression  curves 
may  fit  the  data  with  nearly  the  same  correlation 
coefficient  value,  but  they  give  different  Sg  val¬ 
ues  at  200  cm.  In  any  event,  since  the  regression 
curves  (Fig.  11  and  14)  lie  well  within  each  other's 
standard  deviation,  there  is  no  statistically  sig¬ 
nificant  difference  between  them.  Given  this  rea¬ 
soning,  it  would  seem  that  the  Arctic  and  Antarc¬ 
tic  sea  ice  Sg  vs,  Tp  data  can  be  combined. 

Also  shown  in  Figure  14  is  the  curve  derived 
from  unmodified  eq  1  (Fig.  12).  As  with  the  Beau¬ 
fort  Sea  data,  the  curve  does  not  fit  the  data. 
Taking  the  regression  curve  (Fig.  14)  value  for 


Sgp  =  5.41  %o  at  Tp  =  200  and  dividing  this  salinity 
by  =  34 %o  gives  an  Sr  of  0.159.  Using  this 
value  in  eq  1  gives  the  lower  curve  in  Figure  15, 
Varying  coefficient  a  at  the  Tp  increments  shown 
in  Figure  15  causes  the  Ryvlin  equation  to  give  Sg 
vs.  Tp  values  that  match  the  regression  curve 
through  the  Antarctic  data  (Fig.  14).  Both  the  modi¬ 
fied  Ryvlin  and  regression  curves  are  represented 
by  the  upper  curve  in  Figure  15. 

The  fact  that  growth  rate  coefficient  a  needs 
to  vary  with  Tp  is  not  surprising,  Ryvlin  indicated 
as  much  in  his  description  of  a  in  eq  1.  A  further 
indication  of  the  need  to  vary  a  can  be  seen  in 
Figure  8.  Here  the  time  to  grow  successive  incre¬ 
ments  of  ice  is  shown  to  increase  with  thickness. 
In  short,  the  heat  flux  through  growing  sea  ice 
decreases  exponentially  with  thickness  and  like¬ 
wise  reduces  the  rate  of  ice  growth  (Makshtas 
1984,  Maykut  1986)  and  brine  entrapment.  The 
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Figure  14,  Antarctic  first-year  sea-ice  bulk 
salinity  vs.  floe  thickness  data  of  Jeffries 
(1994).  The  upper  curve  is  the  regression 
fit  to  the  data.  The  lower  curve  represents 
the  results  using  Ryvlin's  equation  for 
ice  grown  in  34%o  seawater. 


Figure  15.  Ryvlin's  equation  for  calculat¬ 
ing  sea-ice  bulk  salinity  vs.  floe  thickness, 
where  parameter  a  is  constant  (lower  curve) 
and  varied  (upper  curve)  to  match  the  re¬ 
gression  curve  through  the  Antarctic  sea- 
ice  data  in  Figure  14. 
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Figure  16.  Sea-ice  stable  salinity  and  the 
sea-ice  seawater  salinity  ratio  vs.  growth 
rate. 


9 


Figure  17.  Antarctic  first-year  sea-ice  bulk  salinity  vs.  floe 
thickness  data  collected  from  several  sources. 


Figure  18.  All  Arctic  and  Antarctic  first-year  sea-ice  bulk 
salinity  vs.  floe  thickness  data  compiled  from  numerous  sources. 


latter,  of  course,  affects  the  salinity  of  the  ice.  This 
is  illustrated  in  Figure  16,  which  shows  the  stable 
salinity  of  sea  ice  vs.  growth  rate.  The  data, 
from  measurements  made  on  Arctic  and  Antarc¬ 
tic  sea  ice  (Nakawo  and  Sinha  1981,  Wakatsuchi 
1983),  clearly  reveal  that  the  salinity  of  sea  ice 
decreases  exponentially  with  decreasing  growth 
rate.  Therefore,  the  growth  rate  coefficient  for  natu¬ 
ral  sea  ice  is  not  a  constant,  but  varies  with  ice 
thickness.  The  scale  on  the  right  side  of  Figure  16 
shows  the  salinity  ratio  of  the  ice  vs.  growth  rate. 
As  expected,  this  ratio  is  not  a  constant,  but  also 
varies  with  ice  salinity. 

In  Figure  14  only  the  Antarctic  sea  ice  data  of 
Jeffries  (1994)  were  used,  to  eliminate  any  varia¬ 


tions  in  sampling  techniques  that  may  exist 
between  different  research  groups  and  to  en¬ 
sure  that  the  ice  sampled  grew  in  seawater  of 
the  same  salinity.  In  the  comparison  between 
this  data  set  and  the  Beaufort  Sea  data,  it  was 
found  that  the  effect  of  seawater  salinity  on 
ice  salinity  was  not  statistically  discernible. 
Therefore,  it  seems  reasonable  to  combine 
the  Antarctic  Sg  vs.  Tp  data  of  Jeffries  (1994) 
with  the  Antarctic  data  of  Fedotov  (1973), 
Wakatsuchi  (1983),  Gow  et  al.  (1987), 
Makshtas  (1984),  Eicken  (1992),  Gow  et  al. 
(1992),  Ackley  et  al.  (1993),  and  Veazey  et  al. 
(1994).  These  data,  315  points,  are  plotted  in 
Figure  17.  There  is  considerable  scatter  in  the 
data,  which  can  be  due  to  a  number  of  fac¬ 
tors:  natural  growth  rate  variations  (Fig.  16), 
brine  drainage  during  coring  and  process¬ 
ing,  errors  in  salinity  determination,  varia¬ 
tions  in  ice  type  (infiltrated  snow,  frazil,  plate¬ 
let  and  columnar  ice),  and  morphology,  such 
as  layered  rafted  ice.  The  salinity  errors  could 
arise  from  errors  in  measuring  the  ice-core 
temperature  and  melt  conductivity,  measure¬ 
ments  that  are  required  for  calculating  ice 
salinity.  The  regression  curve  through  the  data 
in  Figure  17  probably  gives  a  better  assess¬ 
ment  of  Antarctic  sea  ice  Sg  vs.  Tp  than  the 
curve  in  Figure  14  because  of  the  increase  in 
the  number  of  data  points  from  74  to  315. 

Combining  now  the  data  in  Figure  17  with 
the  Arctic  data  in  Figure  11  along  with  first- 
year  sea  ice  S5  vs.  Tp  data  collected  at  Thule, 
Greenland,  and  Hopedale,  Labrador  (Weeks 
and  Anderson  1958),  Fram  Strait  (Gow  et  al. 
1987),  Port  Clarence,  Alaska  (Kovacs  and 
Kalafut  1977),  and  Elson  Lagoon,  Pt.  Barrow, 
Alaska  (L.H.  Shapiro'*'  1994,  personal  contri¬ 
bution,  eight  measurements),  the  distribu¬ 
tion  shown  in  Figure  18  is  obtained.  Interestingly, 
the  regression  curve  through  the  combined  data 
gives  nearly  the  same  Sg  vs.  Tp  values  as  the  curve 
in  Figure  17.  For  example,  at  200  cm  both  equa¬ 
tions  give  5.07%o  for  S5,  while  at  15  cm  the  differ¬ 
ence  is  only  0.30  %o.  This  is  well  within  experi¬ 
mental  error,  and  suggests  that  for  most  purposes 
the  bulk  salinity  vs.  thickness  of  undeformed  first- 
year  sea  ice  grown  in  Arctic  and  Antarctic  seawa¬ 
ter  with  a  salinity  of  about  31  to  34  %o  will  be  very 
similar.  Refinements  and  standardization  of  sam¬ 
pling  techniques  may  one  day  allow  us  to  see  the 


*  Geophysical  Institute,  University  of  Alaska,  Fairbanks, 


10 


Figure  19.  Comparison  between  the  regres¬ 
sion  curve  passed  through  all  the  bulk  salin¬ 
ity  vs.  icefloe  thickness  data  in  Figure  18  and 
the  field  data  in  Figure  7.  The  squares  mark 
the  locations  of  the  ice  thicknesses  used  in  the 
regression  equation  given  in  Figure  18. 


Figure  20.  Beaufort  Sea  multiyear  ice  bulk  salinity  vs.  floe 
thickness. 


effect  of  seawater  salinity  variation  as  well  as 
first-year  sea  ice  morphology  and  structure  dif¬ 
ferences  on  the  bulk  salinity  of  the  ice.  The  litera¬ 
ture  certainly  shows  that  this  can  be  done  under 
controlled  conditions  (Cox  and  Weeks  1975,  Weeks 
and  Lofgren  1967)  and  during  focused  field  mea¬ 
surement  programs  where  the  highly  variable  ther¬ 
modynamic  conditions  controlling  sea  ice  growth 
are  measured  and  accounted  for  (e.g.,  Nakawo 
and  Sinha  1981,  1984;  Wakatsuchi  1983;  Eicken 
1992).  In  the  interim.  Figure  18  provides  a  good 
indication  of  the  bulk  salinity  of  first-year  sea  ice 
vs.  thickness. 

It  may  be  of  interest  now  to  see  how  the  regres¬ 
sion  curve  in  Figure  18  fits  the  data  in  Figure  7. 
Using  the  five  ice  thicknesses  given  in  Figure  7 
and  the  equation  in  Figure  18,  the  bulk  ice  sheet 
salinity  was  calculated.  These  values  are  shown 


as  boxes  in  Figure  19  with  an  interconnecting 
curve.  The  agreement  is  good,  considering  that 
ice  sheet  flooding  in  January  affected  ice 
growth  and  morphology. 

The  results  in  Figure  18  are  representative 
of  first-year  sea  ice  during  the  growth  season. 
In  the  melt  season,  ice  salinity  decreases  due 
to  accelerated  brine  drainage  and  flushing. 
Disregarding  impurity  effects,  the  speed  and 
intensity  of  the  melt  process  are  highly  de¬ 
pendent  on  thermodynamic  factors  such  as 
air  temperature,  solar  radiation  intensity,  snow 
cover,  and  rain,  as  well  as  the  thickness  of  the 
ice.  These  factors  govern  how  fast  meltwater 
forms  on  the  ice  surface  and  the  speed  at  which 
solar  radiation  melts  the  ice  aroxmd  brine  chan¬ 
nels,  pockets,  and  flaws  to  help  open  up  the 
drainage  networks  leading  to  the  sea.  As  indi¬ 
cated  in  Figure  9,  spring  melt  conditions  can 
cause  the  bulk  salinity  of  the  ice  to  decrease  by 
50%  in  one  month.  At  high  latitudes,  where  sum¬ 
mer  melting  is  not  as  intense,  the  process  of  ice 
freshening  is  slower.  Here,  more  than  one  melt 
season  may  be  required  to  achieve  the  reduction 
in  ice  salinity  shown  in  Figure  9.  In  lower  lati¬ 
tudes,  the  melt  season  is  longer  and  sufficiently 
intense  to  completely  melt  the  ice  cover.  Because 
the  degree  of  ice  freshening  in  one  season  is  highly 
variable,  the  salinity  of  second-year  and  older  ice 
also  varies,  but  the  trend  is  decidedly  to  a  lower 
bulk  salinity  for  older  ice. 

Bulk  salinity  vs.  second-  and  multiyear  ice  floe 
thickness  data  collected  in  the  Beaufort  Sea  by 
Kovacs  and  Mellor  (1971),  Cox  (1972'*‘),  Kovacs  et 


*  Only  published  as  shown  in  Figure  6. 
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Figure  21,  Beaufort  Sea  and  Antarctic  multiyear  sea-ice 
bulk  salinity  vs.  floe  thickness. 


Figure  22.  All  Arctic  and  Antarctic  multiyear  sea-ice  bulk 
salinity  vs.  floe  thickness  data. 


al.  (1973),  Cox  and  Weeks  (1974),  Kovacs  (data  on 
file),  and  Meese  (1989)  are  shown  in  Figure  20. 
The  obvious  trend  is  to  lower  S5  values  with  in¬ 
creasing  ice  thickness. 

When  the  higher  S3  vs.  Tp  data  for  Antarctic 
multiyear  ice  (Gow  et  al.  1987)  are  added  to  the 
Beaufort  Sea  data  the  curve  through  the  com¬ 
bined  data  shifts  up  slightly  for  the  thicker  ice 
(Fig.  21).  The  higher  Antarctic  S3  values  may  be 
due  to  ice  structure,  thermodynamic  differences 
during  the  melt  season,  or  differences  in  seawater 
salinity.  When  the  denser  brine  drains  out  of  the 
sea  ice  during  the  melt  season,  some  of  the  drain¬ 
age  system  can  be  reoccupied  by  seawater  (Nied- 
rauer  and  Martin  1979).  In  the  Antarctic,  the  sea¬ 


water  has  a  higher  salinity  than  in  the  southern 
Beaufort  Sea.  Therefore,  the  Antarctic  sea  ice 
could  have  a  higher  salinity  due  to  this  pos¬ 
sible  exchange,  if  all  other  parameters  are  the 
same. 

When  multiyear  sea  ice  S3  vs.  Tp  data 
from  Fram  Strait  (Gow  et  al.  1987)  are  added, 
the  distribution  changes  again,  as  shown  in 
Figure  22.  The  scatter  in  the  data  between  200 
and  300  cm  also  increases.  Some  of  the  lower 
salinity  values  from  Fram  Strait  could  have 
resulted  from  sampling  lower-salinity  melt-pool 
ice.  Whatever  the  reason  for  the  lower  S3  val¬ 
ues,  Figure  22  does  indicate  that  sea  ice  that 
has  survived  one  or  more  melt  seasons  can 
freshen  to  where  the  bulk  salinity  is  1  to  2  %o. 

SUMMARY 

Sea  ice  salinity  data  collected  in  Arctic  and 
Antarctic  waters  indicate  that  the  bulk  ice-sheet 
salinity  decreases  with  increasing  ice  floe  thick¬ 
ness  following  a  negative  exponential-like 
trend.  New  ice  about  0.05  m  thick  has  a  bulk 
salinity  of  about  25  %o.  The  bulk  salinity  of  the 
ice  decreases  rapidly  to  about  6%o  at  an  ice 
thickness  of  about  0.5  m.  After  this,  the  bulk  ice 
salinity  continues  to  decrease,  but  at  a  highly 
reduced  rate,  to  about  5%o  at  the  end  of  the 
growth  season  for  ice  more  than  1.5  m  thick. 
Reasons  for  the  salinity  decrease  are  related  to 
ice  structure,  rate  of  ice  growth,  and  brine  mi¬ 
gration  processes.  The  first  and  second  reasons 
are  related,  in  that  as  the  ice  thickness  increases, 
the  growth  rate  slows  because  the  heat  ex¬ 
change  from  the  ice  bottom  to  the  atmosphere 
is  reduced.  Slower  growth  results  in  the  forma¬ 
tion  of  larger  freshwater  ice  platelets  and  larger 
ice  crystals  with  proportionately  fewer  brine  in¬ 
clusions.  Thus,  the  ice  salinity  is  lower.  Brine  mi¬ 
gration  processes  lead  to  the  gradual  drainage  of 
brine  from  the  ice  to  the  sea  and  a  freshening  of 
the  ice. 

Significant  scatter  and  a  paucity  of  data  pre¬ 
cluded  establishing  a  definitive  relationship  be¬ 
tween  seawater  salinity  and  ice  sheet  structure 
and  the  bulk  salinity  of  the  sea  ice. 

A  significant  decrease  in  ice  floe  bulk  salinity 
occurs  during  the  melt  season.  Depending  on 
many  factors,  including  the  solar  intensity  during 
the  melt  season,  the  bulk  salinity  of  a  first-year 
sea-ice  floe  can  decrease  from  about  5%o  to  about 
2%o.  This  transition  to  second-year  ice  separates 
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the  two  ice  types  according  to  their  bulk  salinity. 
For  this  reason  the  two  ice  types'  bulk  salinity  vs. 
thickness  data  should  not  be  combined,  as  pre¬ 
viously  reported,  in  an  analysis  of  ice  floe  salin¬ 
ity-thickness  trends. 

The  Ryvlin  relationship  for  determining  the 
bulk  salinity  vs.  ice  floe  thickness  is  difficult,  at 
best,  to  use  since  the  growth  rate  parameter  is  not 
a  constant,  but  is  shown  to  vary  with  ice  thick¬ 
ness. 

The  bulk  salinity  vs.  ice  floe  thickness  trend  for 
winter  sea  ice  is  well  represented  by  the  expres¬ 
sion  Sg  =  4.606  +  91.603/Tp  shown  in  Figure  18. 
This  expression  is  applicable  to  both  cold  Arctic 
and  Antarctic  first-year  sea  ice  grown  in  seawater 
of  '-31  to  34  %o  salinity.  The  expression  should  not 
be  used  for  such  southern  sea  ice  environments  as 
the  Gulf  of  St  Lawrence,  where  winter  freeze- 
thaw  cycles  are  common,  and  of  course  the  low- 
salinity  Baltic  Sea,  where  ice  salinities  are  typi¬ 
cally  less  than  1  %o  for  ice  as  thin  as  0.2  m. 
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